Nuclear Charge Radius of ^^Be 
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The nuclear charge radius of ^^Be was precisely determined using the technique of coUinear laser 
spectroscopy on the 2si/2 2pi/2,3/2 transition in the Be^ ion. The mean square charge radius 
increases from ^"Be to ^^Be by 5(r^)"'^^ = 0.69(5) fm^ compared to (5(r-?)"'" = 0.49(5) fm^ for 
the one-neutron halo isotope '^^Be. Calculations in the fermionic molecular dynamics approach show 
a strong sensitivity of the charge radius to the structure of ^'^Be. The experimental charge radius is 
consistent with a breakdown of the N — 8 shell closure. 

PACS numbers: 21.10.Ft, 27.20.+n, 42.62.Fi, 21.60-n, 32.10.Fn 



The explanation of the magic numbers was a first re- 
markable success of the nuclear shell model. In stable 
nuclei the shell structure and therewith associated shell 
gaps are fairly well understood. However, in case of ra- 
dioactive nuclei, far from the /3-stability line, magic num- 
bers can disappear and new magic numbers may arise. 
This was first noticed in case of the = 20 shell closure 
in the region of the neutron-rich Na and Mg isotopes - 
now called the island of inversion. Later it was found 
that this region is not the only island and similar obser- 
vations were made at other locations within the nuclear 
chart (see [ij for a review). The magic number = 8 
is a special case since it naturally arises in every mean- 
field description of nuclei. On the other hand, it is well 
known that light nuclei often exhibit a strongly clustered 
structure and mean-field calculations are not the best 
approach for such nuclei. Hence, a disappearance of the 
N = 8 shell closure might be caused by the cluster struc- 
ture in this region. Moreover, such light nuclei are in 
principle within the reach of ab-initio models using real- 
istic two- and three-body interactions. However, in the 
case of ^^Be only more phenomenological models, e.g., 
the shell model Q, three-body models 0, 0], a micro- 
scopic cluster model a two-center cluster model Q 
and antisymmetrized molecular dynamics are avail- 
able, unfortunately without predictions for the charge 
radius. 

"'^^Be is a key isotope in the beryllium chain, as it is 
located between the one-neutron halo ^-'^Be and the two- 
neutron halo ^^Be. ^^Be has an abnormal ground state 
parity of 1/2+ g, Q indicating the breakdown of the 
N = 8 shell closure with the Si/2 halo orbit lower in 
energy than the pi/2-orbit. Therefore it is expected that 
sd-shell configurations also have considerable influence 
in ^^Be. Such an intruder-state mixing was observed in 



nuclear reaction experiments [iS Eli and a dominating 
sd configuration, i.e., a quenching of the magic shell clo- 
sure at = 8 was postulated. It was argued that the 
lowering of the ^5 /2 level could also cause prolate defor- 
mation (111] and thus a strong clustering effect inside the 
nucleus leading to an extended spatial structure of ^^Be. 
In principle, the intruding s-wave states can also cause 
halolikc structures, but neither the relatively large two- 
neutron separation energy of ^^Be {S2n = 3.67 MeV 12[) 
nor the wide momentum distribution measured in nuclear 
breakup reactions pJi] show the key characteristics of a 
halo structure. However, proton-scattering experiments 
in inverse kinematics revealed a low-density tail in the 
matter distribution, indicating a slight halolike charac- 
ter of the neutron distribution [14]. This is in direct 
contrast with the analysis of ^^Be {d,p) reactions, which 
indicated that the ground state has a small s-wave con- 
tribution (spectroscopic factor of 0.28+007), whereas the 
long-lived excited 0^ state may exhibit an extended den- 
sity tail corresponding to a much larger s-wave spectro- 
scopic factor (O.lStoil) 

[l5| . In this situation additional 
information about the structure of ^^Be is important and 
can be attained by a determination of the nuclear charge 
radius as presented in this paper. 

Charge radii of light nuclei obtained from laser spec- 
troscopy were found to exhibit clear signatures of clus- 
tering effects. This was found in the is otop e chains of 
lithium [3], beryllium [l3] and neon [19|, where the 
charge radii are in good agreement with fermionic molec- 
ular dynamics (FMD) calculations 2^, 21|. In the case of 
beryllium, an extension of our previous calculations for 
^~^^Be towards ^^Be indicates that the nuclear charge 
radius should be very sensitive to the amount of sd-shell 
admixture into this j?-shell nucleus. The charge radius of 
^^Be should be significantly smaller than for ^^Be if the 
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ground state is dominated by a configuration, whereas 
a dominance of an (sd)^ configuration should lead to a 
larger radius. Hence, the observation of an increase in 
the charge radius would directly indicate a quenching of 
the iV = 8 shell. 

We have previously reported on isotope shift 
measurements for the beryllium isotopes ^~^^Be at 
ISOLDE/CERN [11, [2^. A measurement on ^^Be turned 
out to be unfeasible at that time because of its low pro- 
duction rate. In this letter, we present such a measure- 
ment using an improved and more sensitive technique 
and we compare the extracted charge radius with new 
FMD calculations. 

We measured the optical transition frequencies of ^^Be 
ions for the excitation 2sij2 — > '^Pi/2,3,/2 (Dl and D2 
lines) at 313 nm using the technique of frequency- comb 
based collinear and anticollinear laser spectroscopy on 
a Be+ beam, which is independent of inaccuracies in 
calibrating the ion beam velocity [22|. Details of the 
experimental setup were presented in fl7| . Hence, we 
give here a brief description, highlighting particularly the 
changes made for the ^^Be measurement: ^^Be nuclei 
were produced at ISOLDE/CERN by pulses of about 
3 • 10^'^ protons at 1.4 GeV impinging on a UC^ tar- 
get. On average about 1 pulse every 3-4 s was hitting 
the target which is coupled to the RILIS resonance ion- 
ization laser-ion source. Typically 8000 -'^^Be ions were 
extracted within 60 ms after the proton pulse, which is 
three times the half-life. The 40 keV ion beam was su- 
perimposed with copropagating and counter-propagating 
laser beams. Doppler-tuning was performed by applying 
a post-acceleration voltage to the fluorescence detection 
region (FDR). Behind the FDR, the ions were deflected 
and detected with a secondary electron multiplier tube 
(SEM). To suppress background from laser stray light, 
an ion-photon coincidence technique Q was used. For 
this purpose the signals from the photomultiplier tubes 
were digitally delayed using a field programmable gate 
array (FPGA) and only pulses from the photomultiplier 
in coincidence with the corresponding ion signal were 
recorded. This avoids dead time losses from the elec- 
tronic delay circuit that compensates for the time of flight 
of the ion from the FDR to the SEM. Only pulses from 
the photomultiplier in coincidence with the correspond- 
ing ion signal were recorded. 

Typically, the resonance spectra in the Dl and D2 tran- 
sitions of -'^^Be+ ions were recorded with a signal to noise 
{S/N) ratio of about 10 after about 2 hours. In total 
25 and 15 spectra for each laser beam direction were 
recorded for the Dl and D2 transition, respectively. Since 
^^Be has no hyperfine structure these spectra were fitted 
with a single Voigt profile as shown in Fig. [TJ The volt- 
age scale was converted to the respective frequency scale. 
The full width at half maximum of a resonance is about 
40 MHz, composed of a Lorentzian contribution of typ- 
ically 25 MHz, which is close to the natural linewidth 
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FIG. 1: (Color online) Fluorescence spectra of ^^Be in the 
Dl transition as a function of the Doppler-tuned frequency 
in collinear (left) and anticollinear (right) direction. A Voigt 
profile is fitted to the data points. The lower trace shows the 
residua of the fitting. 

of about 20 MHz, and a Doppler width of 30 MHz. No 
asymmetry in the line shape of any isotope was observed 
in contrast to earlier measurements [2^. With the peak 
positions determined for both laser directions {va, Vc) the 
transition rest frame frequency was calculated accord- 
ing to 

vl = Va- Vc- (1) 

The absolute frequencies are then used to determine 
the isotope shifts 5v^^ = i'o(^Be) — i'o(^Be). Compared 
to our previous work, we could considerably reduce the 
main source of systematic uncertainty caused by a pos- 
sible photon recoil shift: The even isotopes ^"Be and 
^^Be have no hyperfine structure and provide a closed 
two-level system which can repeatedly scatter photons 
of the laser beam. Resonant interaction with the anti- 
collinear laser beam decelerates the ions while the inter- 
action with photons from the collinear beam accelerates 
them. This shifts both frequencies, Va and Uc, to higher 
values and does not cancel when applying Eq. ([T]). We 
investigated this effect using ^°Be and determined the 
absolute frequency for laser powers varying from 0.5 to 
8 mW with a laser beam diameter of about 3 mm. No 
systematic shift as a function of the laser power was ob- 
served within a standard deviation of 200 kHz. With a 
typical laser power of 3 mW, an upper limit of 100 kHz 
is assumed here as a systematic uncertainty. Addition- 
ally, a possible laser-ion beam misalignment contributes 
with an upper limit of 500 kHz according to our pre- 
vious experiments. The isotopes ^-I'^'i^Be were reinvesti- 
gated in the Dl line and ^°Be additionally in the D2 line. 
The new measurements on ^^Be were interspersed with 
measurements on ^"Be throughout the beamtime, in or- 
der to check them with a well-known reference. Isotope 
shifts 5v^^^ ^ as listed in Tab. H were calculated using 
the absolute frequency of ^Be as obtained as the aver- 
age of all beamtimes for the respective transition. The 
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FIG. 2: (Color online) Experimental charge radii of beryllium 
isotopes from isotope shift measurements in 2008 (o) and in 
2010 (this study •) compared with the FMD results (this work 
■ ). Charge radii calculated for the pure p^-configuration (V) 
and pure (sd)^ intruder configuration (A) are also indicated. 



isotope shifts determined in 2008 and 2010 agree within 
300 kHz. The change in the mean square charge radius 
is calculated according to 
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with the calculated mass-shift values (5j^ms ' partially up- 
dated for more accurate mass measurements [2^, [2^ , and 
the electronic field-shift factor F = -17.02 MHz/fm^, 
both from [H, [IBj. While the mass-shift of "Be in- 
cludes a 208(21) kHz contribution from nuclear polariz- 
ability [2^, this effect can not be calculated for ^^Be for 
which the low-lying dipole strength distribution is not 
yet known. However, from preliminary breakup reaction 
data it is expected to be smaller than for ^^Be. Hence, 
an additional uncertainty of 60 kHz was added to the 
mass-shift of ^^Be. 

To determine absolute rms charge radii we have used 
the reference radius Rd^Be) = 2.519(12) fm 27| from 
elastic electron scattering. Including the uncertainty of 
the reference radius the charge radius of ^^Be was con- 
sistently extracted, from the Dl and D2 transition to be 
Rci^^Be) = 2.502(16) fm with a relative uncertainty of 
less than 1 %, dominated by the uncertainty of the refer- 
ence radius. The nuclear charge radii are listed in Table U 
and plotted in Fig. [2] Results of our present (•) and our 
previous work (o) are compared with new calculations in 
the FMD approach. 

To understand the consequences of the result for the 
structure of ^^Be, we reinvestigated our FMD calcula- 
tions for the beryllium isotopic chain. FMD [s^l uses a 
Gaussian wave-packet basis for the single-particle states 
that allows to describe nuclei with clustering and halos 



in a consistent picture. Good quantum numbers are ob- 
tained by projecting the intrinsic many-body basis states 
on parity, angular and total linear momentum. The pa- 
rameters of the basis states are obtained by variation 
after projection for the spins of the low-lying states for 
each isotope. Additional basis states are obtained by con- 
straining the quadrupole moment of the intrinsic states. 
The full wave functions are obtained in a multiconfigu- 
ration mixing calculation with 30-40 basis states. 

In [23| we performed FMD calculations that were able 
to describe the evolution of the charge radii in the beryl- 
lium isotopes up to ^^Be. However with the effective 
interaction employed there, the gap between p- and sd- 
shell-orbits is too large and the parity inversion in ^^Be 
could not be described. As the competition between p- 
and sc?-configurations plays an essential role in ^^Be, a 
different interaction is used in this work. It is derived 



from the Argonne V18 potential [3l| in the unitary cor- 
relation operator method 3^ with correlation functions 
obtained from a Hamiltonian evolved in the similarity 
renormalization group [s^. This UCOM(SRG) interac- 
tion has been used successfully in FMD for a calculation 
of the ^He(a,7)''Be capture reaction [33|. However no- 
core shell model calculations [s^ show that the effective 
spin-orbit splittings with this realistic two-body interac- 
tion are too small. To correct for that deficiency without 
employing three-body forces we introduce a phenomeno- 
logical parameter rj to change the strength of the spin- 
orbit force in the S* = 1, T = 1 channel. A factor 77 « 2 
works remarkably well for the beryllium isotopes. This 
modification of the spin-orbit force leads to increased 
binding and it shifts the relative position of p- and sd- 
orbits. The best agreement with experimental separation 
energies is obtained for a strength factor which changes 
from 77 = 1.9 (for ^Be) to ry = 2.2 (for ^^Be). With the 
exception of ^^Be properties like radii and transitions de- 
pend only very weakly on the details of the spin-orbit 
force. The A-dependence of 77 indicates that this proce- 
dure might indeed simulate the role of three-body forces. 

The intrinsic states of the beryllium isotopes show a 
more or less pronounced a-cluster structure as already 
discussed in [20|. Depending on the isotope the distance 
between the a-clusters is different. This effect is the rea- 
son for the large decrease in the charge radius from ^Be 
to ^°Be. The jump between ^°Be and ^^Be on the other 
hand, is mainly explained by the motion of the ^°Be core 
with respect to the total center-of-mass of the core and 
the halo neutron in an extended s-orbit. However, we 
also find a small contribution related to core excitation 
from the c?-wave component. In the case of the O"*" states 
of ^^Be we find in the FMD calculation two different lo- 
cal minima that correspond essentially to a p^ and an 
intruder (sd)^ configuration as shown in Fig. [3l Whereas 
the charge radius of the p^ configuration with 2.33 fm is 
only slightly larger than in ^°Be, the (sd)^ configuration 
has a significantly larger charge radius of 2.44 fm, which 
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TABLE I: Changes of mean square charge radii S (r^) calculated from measured isotope shifts Svf^''^ relative to ^Be and 
theoretical mass-shift values Sv^g [2342^ according to Eq. ([2|. J?c is the corresponding rms charge radius based on the 
reference radius of ^Be [IJ]- For comparison, the matter radii Rm [l^l and the results of FMD calculations are given as well. * 
taken from [itI ]. 

Isotope (^t/fs^ (MHz) (5t/^s (MHz) g(rg)(fm^) _Rc (fm) _Re (FMD) (fm) i^^ (Exp.) (fm) (FMD) (fm) 



'Be 


-49 236.88(97)* 


-49 225.779(38) 


0.66(6) 


2.646(16) 


2.59 


2.31(2) 


2.37 


«Be 











2.519(12) 


2.50 


2.38(1) 


2.44 


i"Be 


17 323.45(83) 


17 310.441(12) 


-0.77(5) 


2.361(17) 


2.31 


2.30(2) 


2.23 


D2 if'Be 


17 325.6(16) 


17 312.569(13) 


-0.77(9) 


2.361(24) 








"Be 


31 564.71(73) 


31 560.294(24) 


-0.26(4) 


2.466(15) 


2.38 


2.73(5) 


2.80 


i^Be 


43 391.46(80) 


43 390.168(39) 


-0.08(5) 


2.503(15) 


2.41 


2.59(6) 


2.52 


D2 i^Be 


43 397.0(16) 


43 395.499(39) 


-0.09(9) 


2.501(22) 








"Be 










2.42 


3.16(38) 


2.74 



is caused by both an increased distance of the a-clusters 
and a correlation with the neutrons which leads to an 
a-^He like structure. The (sd)^ configuration in ^^Be 
is very different from the picture in ^^Be where a very 
small neutron separation energy favors the appearance 
of an s-wave halo. In ^^Be the two neutrons are much 
better bound and the FMD single-particle orbits are lin- 
ear combinations of d- and s-orbits. This is also reflected 
in the matter radius that is significantly smaller in the 
^^Be (sd)^ configuration than in ^^Be. As the and 
(sd)^ configurations are very close in energy, the mixing 
is strong and the charge radius of the ground state de- 
pends sensitively on the mixing. This is also the case 
for other observables, like the monopole matrix element 
between the two 0^ states and the quadrupole transition 
strengths of the two 0+ states to the first 2+ state. With 
the strength of spin-orbit force adjusted to the separation 
energy we obtain a good agreement for all transitions as 
shown in Table HIl The (sd)^ configuration contributes 
about 70% to the ground state. This is reflected in the 
charge radius which is larger than the calculated ^^Be 
charge radius by about 0.03 fm. The FMD calculations 
are therefore able to describe the experimental trend of 
the charge radii from ^Be to ^^Be. However, the absolute 
values are somewhat too small, especially for the heavier 
isotopes. This might be related to deficiencies of the ef- 
fective two-body interaction which is known to saturate 
at too high densities for heavier nuclei (ssf . 

While the charge radius increases from ^^Be to ^^Be, 



TABLE H: Two-neutron separation energy and M{EO) and 
B{E2) transition strengths in ^^Be. 







FMD 


Exp 


5'2„(MeV) 
B{E2;2+ - 
B{E2-0+ - 
M{EO- 0+ - 


> 2+)(e2fm4) 
->02+)(efm2) 


3.81 
8.75 
7.45 
0.90 


3.673 \12] 
8.0(3.0) [35] 
7.0(6) [36] 
0.87(3) [36] 




y (fm) y (fm) 

FIG. 3: Cuts through the intrinsic densities of the ^^Be states 
inp^ (left) and (sd)^ configurations (right) obtained by varia- 
tion after projection. The contour lines are in units of nuclear 
matter density. 



the matter radius decreases considerably [28], which is 
also in accordance with our FMD calculations. 

In summary, we investigated the 2si/2 '2pi/2.3/2 
transitions in ^^Be ions using combined coUinear and an- 
ticollinear laser spectroscopy to extract the isotope shifts 
of ^^Be relative to ^Be. By combining the experimen- 
tal data with recent mass-shift calculations we obtained 
the change in the mean square charge radius from ^"Be 
to ^^Be of 0.69(5)fm^ resulting in an absolute radius of 
2.502(16)fm. This increase of the charge radius provides 
important new information for understanding the struc- 
ture of ^^Be. In FMD model calculations the increase of 
the charge radius is related to a breakdown of the = 8 
shell closure with an (sd)^ admixture of about 70%. This 
is also consistent with other observables such as the life- 
time of the isomeric O"*" state and the quadrupole transi- 
tion strengths to the first 2+ state. 
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